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ABSTRACT 


A study and experimental investigation has been performed to determine the 
feasibility of measuring regional blood flow and volume in man by means of micro- 
wave radiometry. An indication was expected of regional blood flow from measure- 
ment of surface and subsurface tenq>eratures with a sensitive radiometer. 

Following theoretical modeling of biological tissue, to determine the 
optlnnin operating frequency for adequate sensing depth, a sensitive microwave 
radiometer was designed for operation at 795 MHz. A temperature sensitivity of 
of 0.06°K rms was realized In this equipment. 

Measurements performed on phantom tissue models, consisting of beef fat and 
lean beefsteak showed that the radiometer was capable of sensing temperatures 
from a depth between 3.8 and 5.1 cm. 

Radiometric and thermodynamic temperature measurements were also performed 
on the hind thighs of large dogs. These showed that the radiometer could sense 
subsurface temperatures from a depth of, at least, 1.3 cm. Delays caused by 
externally-generated RF Interference, coupled with the lack of reliable blood flow 
measurement equipment, prevented correlation of radiometer readings with regional 
blood flow. For the same reasons. It was not possible to extend the radiometric 
observations to human subjects. 

A set of conclusions and recommendations are presented at the end of the 
report. 
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Section 1 


INTRODUCTION 

The purpose of the research described in this report was to convi ict i cudy 
and experimental investigation aimed at determining the feasibility cf measuring 
regional blood flow and volumes in man by means of microwave radiometric sensing 
techniques* Since microwave radiometers are temperature sensirg instruments, it 
was expected that an indication of regional blood flow would be gained by measur- 
ing related temperatures down to a depth dictated by the selected operating 
frequency. The study has confirmed the ability of a sensitive radiometer to sense 
subsurface temperatures in living tissue, although the work schedule did not allow 
sufficient time for correlation of this information with regional blood flo\^7. 

The project involved the following basic elements of work; 

1) A theoretical study to determine radiometric sensing depths of micrc'^;a^c 
thermal energy, emitted by models of biological tissue, over the frequency range 
0.13 to 5 GHz. The purpose of this work was to locate an optimum operating 
frequency for the radiometer. This was established at 725 MHz. 

2) Design and development of a sensitive microwave radiometer suitable for 
laboratory-*t 3 rpe observations on tissue models, animals and human subjects. This 
was successfully accomplished; the radiometer operates at a frequency of 793 MHz 
and has a minimum detectable temperature sensitivity of 0.06 degree Kelvin rms, 
for a 10-second output time constant. 

3) Radiometric and correlative temperature measurements on phantom mouels of 
biological tissue, consisting of beef fat and lean beefsteak, to determine the 
depth from which temperatures can be sensed. It was found that the radiometer 
sensed temperatures from a depth between 3.8 and 5.1 cm in lean beefsteak. 

4) Temperature profile versus blood flow measurements on hind thighs of large 
dogs I without radiometric observations. Correlation of blood flow with surface 
and subsurface temperatures was accomplished with somewhat limited success due to 
problems encountered with the stability of electromagnetic flowmeters. Success 
was achieved only on a "go no-go" basis, in that measurements were made under full- 
flow or no-flow conditions - the no-flow condition was achieved through the use of 


a tourniquet. Some temperature profiles were obtained after a given area wad 
coviled with an ice-alcohol mixture. 

5) Radiometric and correlative temperature observations on the hind thigli 
a large dog. This investigation showed that the radiometer sensed temperat>uci> 
liom a depth of, at least, 1.3 cm. 

The above work demonstrated clearly that radiometric measurement of subsur- 
face temperature changes, associated with changes in blood flow, is feasible. 
Unfortunately, due to time limitations Imposed by design changes, caused by 
externally generated RF interference, and the lack of adequate blood flow measur- 
ing equipment, it was not possible to establish a quantitative relationship 
between radiometer apparent temperatures and regional blood flow. The shortage 
of time also precluded any observations on human subjects. 

Considering the problems encountered with externally-generated RF interference 
and the lack of reliable blood flow measurement equipment, it is felt that a good 
beginning has been made in the application of microwave radiometry to subsurface 
temperature sensing in living tissue* The results of the study are presented in 
the form of numerous plots, tables and related discussions. 

A set of conclusions and recommendations are presented at the end of the 
leport. It is felt lhat, if the recommendations are adopted and the requisite 
observations made on human subjects, the laicrowave radiometric technique can be 
properly developed and ultimately become quite useful for subsurface temperature 
sensing and indirect measurement of blood distribution in the human body. 



Section 2 


THEORETICAL STUDIES 


Preliminary theoretical studies were conducted to determine the feasibility 
of using a radiometer to detect internal blood flew in living tissue. These 
studies also provided a method for selecting the optimum radiometer frequency for 
detecting blood flow. The tissue structure was modeled by a planar layered medium, 
of infinite lateral extent, because such a model allows for direct computations of 
brightness temperature. The temperature of the layer representing the region of 
internal blood fl'^w was assumed to be indicative of the volume of blood flow. The 
optimum radiometer operating frequency for detecting changes in blood flow was 
determined by changing the temperature of this layer and finding the frequency at 
which the largest change in surface brightness temperature occurred. 

2.1 DIEIECTRIC PERMITTIVITIES OF BIOLOGICAL TISSUES 

The tissue structure was modeled by a planar layered medium of infinite 
lateral extent. The five media, shown in Figure 2-1, used for the model were, in 
descending order: the antenna, dielectric, skin, fat, muscle and bone. The antenna 

dielectric and bone layers were semi-infinite i'^' thickness. The thickness of the 
skin layer was 0.3 cm, the fat layer was 0.5 cm, and the muscle layer thickness was 
a variable in the model. Before surface brightness temperatures could be computed, 
the dielectric permittivities of each layer needed to be specified. 

Samples of lean beefsteak, representing muscle, beef bone, abdominal human fat, 
and human skin were submitted to the Laboratory for Insulation Research, at the 
Massachusetts Institute of Technology, for dielectric permittivity measurements at 
0.5, 1 and 3 GHz. The beef bone appeared to be taken near a joint and contained 
little oiarrow. Measurements were made at two temperatures; 56°C and 40*^C for the 
beefsteak, bone, and skin, and 27.5°C and 42.9'^C for the fat. The permittivity of 
fat, at 36‘^C, was determined by a linear interpolation between 27.3'^C and 42.9<>C. 
Additional infemration was obtained from Pennock^ on the permittivity of human skin, 
fat, and muscle a. 0.15 GHz and at a temperature of 37^C. The measurements, made at 
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Figure 2-1 - Phantom Tissue Model 
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two temperatures, indicate that variations in permittivity with temperature is 
slight over a range of 1°C, and Pennock's data at 37°C should not be signif icar*-T v 
different from measurements made at 360C. The permittivity of bone at 0.15 GHz 
was determined by extrapolating the higher frequency data, noting that the varia- 
tion of bone permittivity with frequency is slight. 

Figures 2-2 through 2-5 show the complex permittivity, c' - jt”, as a function 
of frequency for a temperature of 36°C. By interpolating between the four data 
points, the permittivity at any frequency between 0.15 and 3 GHz can be approximate . 
The error in all measurements was about 5%. One factor which could contribute to 
the permittivity of dead tissue being different from living tissue is the absence 
of circulating fluids. Accordingly, great care was exercised in sealing tissue 
samples against evaporation. 

The antenna dielectric was chosen to match, as close as possible, the tissue 
structure, while having negligible attenuation. To match the impedance of the 
model, consisting mainly of 'luscle tissue, in which the real part of the permit- 
tivity is high, the antenna should also have a large real permittivity. In order 
to exhibit negligible attenuation, the antenna imaginary permittivity should be as 
small as possible. The antenna chosen to meet these requirements has a permittivity 
of 30- jo. 015. Further details on the antenna are given in Section 3. 

A parametric study was conducted during the report period t- assess the 
influence of the antenna dielectric constant on tissue brightness te iperatures. The 
parameter involved in the study was fat thickness lying on top of muscle tissue. 

Figure 2-6 shows the results of calculations performed, using Program lAYER, 
at 725 MHz. A muscle tissue thickness of 7 cm was used in the model. The plots 
show that, for a lean individual with a fat layer of 0.5 cm, the optimum dielectric 
constant is approximately 25. As the fat layer increases to 1.5 and 3.0 cm, the 
loading material should exhibir dielectric constants of 7 and 2.5, respectively. 

Since the peaks in the curves are fairly flat, considerable latitude should be 
acceptable in the value of the dielectric constant; this is estiinate at - 40%, for 
a brightness temperature error of 5‘’K. 
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Figure 2-2 - Dielectric Permittivity of Human Skin versus Frequency 
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Figure 2-4 - Dielectric Permittivity of Lean Beefsteak veraua Frequency 
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Figure 2-5 - Dielectric Permittivity of Beef Bone versus Frequency 
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As mentioned above, the antenna furnished for the radiometer is loaded with 
a material having a dielectric constant of 30. Based on the above criterion, this 
antei na should furnish adequate absolute brightness temperatures on all individuals 
exhibiting a fat layer of up to approximately 1.5 cm. This is expected to encompass 
a fairly large percentage of the male pupulation. 

2.2 BRIGHTNESS TEMPERATURE AND EMISSIVITY OF LAYERED MEDIA 

The brightness temperature, Tg, referred to herein, is a measure of power 
being thermally radiated from the tissue layers upward into the antenna. The antenna 
also radiates power downward. A measure of this power is the internal radiometer 
temperature, T^. A portion of Tg is reflected by the tissue layers back into the 
antenna. The apparent temperature, T^, is a measure of the total power travelling 
upward into the antenna and is given by, 

T^ a Tg ♦ RTg , Or (2-1) 

where, R is the power reflectivity of the tissue layers. R is related to the 
emissivity t of the tissue-antenna structure simply by, 


R » 1-e (2-2) 

The emissivity t can be found by assigning a constant temperature, T^, to the layered 
medium, computing the brightness temperature, Tg, and then dividing Tg by T^. That 
is, 

e » Tg/To (2-3) 

Any Tq can be used because Tg is directly proportional to Tq i.e., the relationship 
between Tb and Tg is linear. 

The brightness temperature of a layered medium depends on the permittivities, 
thicknesses, and temperatures of the various layers and on th*:- radiation frequency 
and angle of incidence. With these inputs, a computer program LAYER (developed 
during previously condU(;ted work under Company sponsorship) can compute Tg. 
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Computations at five frequencies (0.15, 0.5, 01725, 1.0, and 3.0 GHz) were 
peituimed fni two muscle layer thicknesses ot 5 and 7 cm. The temperatures of the 
layer interfaces weie, in descending o*.der, 306.0, 306.5, 307.0, 309. 0°K, as shown 
in Figure 2-1. The temperature gradient within the layers was taken to be linear. 
Zero incidence angle was always used. Figure 2-7 shows the results of ttese computa- 
tions. At the higher frequencies, radiation originating deep within the muscle 
layer and from the bone layer is, to a large degree, attenuated by the time it 
reaches the antenna. This explains the convergence of the curves for the 5 cm and 
7 cm models at higher frequencies. The rise in brightness temperature, with increas- 
ing frequency, is due to a better impedance match between the antenna and the remain- 
ing layers at higher frequencies. 

2.5 SELECTION OF RADIOMETER OPERATING FREQUENCY 

In selecting the radiometer operating frequency, the temperature of the muscle 
layer was assumed to be indicative of the blood flow to be detected by the radiometer. 
By changing the muscle layer temperature and noting at which frequencies the largest 
change in brightness temperature occurs, the optimum radiometer frequency for sensing 
blood flow changes can be determined. Computations similar to those described in 
Section 2.2 were conducted in which the temperature of the muscle-bone interface 
was varied in increments of 0.2°K, from 307°K to 311°K. The brightness tempera- 
ture, Tg, varies linearly with the temperature of the bottom interface. The 
change, ATg, in brightness temperature per change, in internal tenqperature 

of the bottom interface is plotted against frequency in Figure 2-8. A maximum 
is exhibited near 0.5 GHz, with the curves being fairly flat over the region of 
0.15 to 1.0 GHz. The power attenuation through the various media varies directly 
with the imaginary part, e", of the permittivity and with frequency. Figures 2-2 
ttuough 2-5 show that, at the lower frequencies, e” displays a sharp rise which 
causes the curves in Figure 2-8 to fall off with decreasing frequency. At the 
higher frequencies, where e" is fairly constant, the ATg/ATo curves also fall off 
because of the dependence of power attenuation on frequency. 

Since the ratio, < Tg/ATf,, is a maximum at a frequency of 0.5 OTz, chis would 
be an obvious choice of ladiometer operating frequency. However, it was felt that 
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Figure 2-8 - Ratio of Change in Brightness Temperature to Change in 
Temperature versus Frequency, for Muscle Tissue Model 
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operation at such a low frequency would result in a rather large antenna aperture 
and, hence, poor spatial resolution. Acccordingly, a search was conducted for a 
clear region higher in the UHF spectrum, referring to the Federal Conmiunications 
Counission frequency allocations. The result of this work showed that the region 
between 650 and 800 MHz was clear, except for the Channel 56 OHF TV channel at 723 
to 728 MHz. It was, therefore, decided to set the radiometer operating >'>and in this 
frequency interval and provide a narrow-band rejection filter to eliminate any 
possible interference from this television channel. 
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Section 3 


RADIOMEIXR DBSl(a« 


3.1 DESIOl OBJECTIVES 

The design objectives for the microwave radiometer were as follovrs; 

1. Operation at a fixed frequency, located somewhere In the range 0.5 
to 5 GHz. 

2 . High temperature sensitivity - 0.06 degree Kelvin for a ten (10) 

second output time constant (or 0.2OK for a one (1) second time constant), 

3. A manually controlled, movable RF head with a specially designed, small 
waveguide horn antenna, and 

4. A portable console containing the low frequency (video) portion of the 
radiometer, and a signal processing and recording unit. 

The above objectives have been realized with a relatively lightweight and 
portable system suitable for laboratory Investigations of blood flow. 

5.2 BASIC RADIOMETER 00NF16URATI(»( 

2 

The basic configuration of this equipment Is that of a load-comparison Dlcke 
radiometer. This configuration greatly reduces spurious noise signals, caused by 
short-term receiver gain fluctuations, by shifting the post-detection sl^pial spectrum 
up from DC to a carrier which Is above the highest expected gain variation frequency. 
This is accomplished by modulating the RF signal, at a point between the antenna and 
receiver, and coherently detecting this signal after the second detectcr. It can be 
shown that changes In radiometer output voltage, due to short-term variations in 
receiver gain, are zero. 

The modulating waveform la a 50%-duty-cycle square-wave; thus, the radiometer 
dwells 50% of the time on the antenna signal and 50% of the time on a resistive 
termination, commonly referred to as the reference load. This comparison between 
the input signal and reference load is performed by a diode-type switch. 
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Figure 3-1 shows a photograph of the radiometer. The individual components 
are shown labelled in Figure 3-2. A functional block diagram of the radiometric 
system is presented in Figure 3-3. Basically, the system consistr of an antenna; 
internal calibration section; modulator switch (Dicke load-compa: ^un switch); 

RF amplifier section; square-law detector/video amplifier unit; lock-in amplifier 
with coherent detector and low pass output filters; and a signal processing and 
record5.ng unit. Further details, concerning these subdivisions, are presented in 
Paragraphs 3.3 and 3.4. 

3.3 RADIOMETER DES1(^ AND PERFORMANCE CHARACTERISTICS 

The design and performance characteristics of the radiometer are presented 
in the paragraphs to follow in terms of minimum detectable temperature sensitivity 
and other important data. 

3.3.1 RADIOMETER TEMPERATURE SENSITIVITY 


The minimum detectable temperature sensitivity of any radiometer is usually 
expressed in terms of a unit signal-to-noise power ratio at the output of the 
radiometer. This definition is based upon a number of parameters whose actual 
values are referred to the input of the radiometer. These parameters include the 
input losses, input noise figure, and the effective passband of the radiometer. 


The minimum detectable temperature sensitivity, AT, is usually expressed 
as follows: 


{(PL-I) 




, rms 


(5-1) 


where, 

k is a constant depending on the character of the waveform employed for 
input signal commutation; this is usually set at a value of 2 for 
square-wave input chopping. 

(FL-1) represents the effective receiver noise figure. 
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Figure 3-2 - Major Consponenta in UHF Mlcrewave Radiometer 
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Figure 3-3 - Functional Block Diagram of UHF Radiometer 
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F is the noise figure of the first detector (in this particular radioaeter 
it is the first RF amplifier). 

L represents the sun of all the input losses, between the mouth of '.^e 
antenna and the first RF anplifier. 

is the ambient radiometer RF head temperature (*^K). 

is the average apparent temperature seen by the antenna (OK^ . 

B is the 3-db RF bandwidth prior to detection (Hz). 

T is the radiometer output time constant (seconds). 

In this design, the above parameters are as follows: 
k s 2 

F s 2.4 db (nuHterical value » 1.74) 

L * 2.2 db (numerical value ■ 1.66) 

T ■ 3110K 

o 

T « 200°K 

a 

B « 75 MHz at NASA--ADie8 • B » 119 MHz at Radioaetric Technology » lnc« 

Fron the above values it is possible to calculate the tesqperature sensitivity, 
AT, for an output time constant of one (1) second: 


. 2ra.7*X 1.66 - » 311 * 200] . „ ^ 

^5 X 10® X 1 

For a 10-second tisie ccHistant, the temperature sensitivity would iiiq>rove by 

a factor of Jio , or 3.16. Thus, for this time constant, _ 0.19 _ ^ 

ni ■ ^ ■ ■ o.oo^ ms 

3.16 

These sensitivities meet the original objective of 0.2°K and 0.06®K, for time con- 
stants of 1 and 10 seconds, respectively. 
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3.5.2 OTHER DBSICa; AIQ) PERFORMAHC£ CHARACTERISTICS 

The other important radiometer performance characteristics are summarized 

below: 

Operating RF center frequency 793 MHz at NASA-Ames 

V65 MHz at Radiometric Technology, Inc. 

Antenna type Dielectric-loaded waveguide horn antenna. 

Dielectric constant of antenna 30-j0.015 (Emerson A Cunminga, Inc. 

loading material stycast). 

Antenna over-all dimensions 2.220 X 1..142 in. 

Antenna inside dimensions 2.156 x 1.076 in. 

Teflon antenna boot 1/32 in. teflon, dielectric constant a 

2.1 - jO.017. 

Antenna and input losses 2.2 db total. 

Type ot RF amplifier Low noise transistor amplifier, operating 

at ambient temperature. First anq>lifier 
noise figure = 2.4 db. 

RF amplifier gain (2 amplifiers) 23 43 a 66 db. 

RF 3-db bandwidth 75 MHz at NASA-Ames; 119 MHz at Radio- 

metric Technology, Inc. RF bandwidth 
defined by bandpass filter. 

Modulator switch switching wave 1,000 Hz 50%-duty-cycle square wave. 

Output detector/anq>lif ier Square-law characteristic in detector. 

Video amplifier integral with detector. 
Amplifier gain ■ 20 db. 

Input calibration method 2-point calibration, at 0 and 53°R, 

using noise diode source. 

RF head temperature Stabilized at 311 i 0.5OR after 45- 

minute warm-up. 

Output time constants 0.1, 0.3, 1, 3, 10 A 30 secs., and 

external connection. Also, 1, 3, 10 A 30 
msecs, and external connection. These 
are available in the Lock-In Amplifier unit. 
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Badloneter temperature aensitlvity 

Radiometer absolute accuracy 

Signal processing and recording 
method 

Radiometer input power 
Auxiliary instrumentation 


3.4 RADIOMEISR DESIGN DEIIAIIS 


0.18°K rti'.s, for 1 sec., time cons cant 
0.06°K rms, for 10 sec., ..imi <'onstauc. 

± 0.5®K, at or near a null ft the input. 

Analog-to-digital conve'-sim of radio- 
meter and tenq>err"ure pi jbe inputs. 
Recorded on 2-track incremental digital 
magnetic tape cassette recorder. 

115 VAC, 60 Hz, 2.5 amps., including 
heater. 

Antenna calibration box. Tl.ermistor 
temperature probes; 

Surface temperature 
Meedle-tvpe temperature probe^ 

Rectal temperature 
Air temperature. 


The design aspects of the various subdivisions of the radiometer will be 
discussed with reference to Figure 3-3, which shows the functional block diagram 
of the system. 

3.4.1 DIELECTRIC-LOADED WAVEGUIDE ANTENNA 


A factor that affects the transfer of power from the human body to an open 
waveguide antenna, connected to a radiometer, is the rather large disparity between 
the characteristic impedance of the body and that of free space. If one were to 
use a horn which is designed to produce efficient transfer of power from free space 
into a waveguide, as a means for transferring power from the huaan body to the wave- 
guide, one would be confronted with a somewhat large impedance mismatch. 

The characteristic impedance Zqq of free space is equal to ■ 37 / ohms. 

•^0 

The characteristic impedance of muscle tissue, at 800 MHz can be determii cd from 
the fact that its dielectric constant, K, is approximately 58- j38. Substitu'-lng 
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this value of K into the formula for the characteristic impecance, It is found 
that it is approximately one-seventh that of free space and that, because of the 
finite conductivity, this impedance is not quite real. A more important factor 
is that the impedance, Z^, of the body is relatively low compared with that of 
free space. 

If one were to design p horn, for transfer of power from a waveguide into 
free space, and vice versa, and fill it with a material ' . lag a dielectric constant 
of approxiaiately 58 and, at the same time, iacrease the dimensions of the horn by 
the factor JSS, the horn would transfer power efficiently into a medium with a 
dielectric constant of 58. 

The dimensions of the waveguide horn would shrink by a factor of about seven. 
Such a horn can be fed with a coaxial line. Since the characteristic impedance of 
the waveguide, filled with a dielectric K *■ 58, would be about one-seventh of the 
conventional waveguide, the interface mismatch between the 50-ohm coaxial line and 
this device would be less than it is usually; thus the device for coupling a coaxial 
line to this type of an anthnna does not present a design problem. 

The antenna designed for the radiometer reflects the above approach. Figure 
3-4 shows two views of the antenna. It will be noted that the dielectric material, 
used for loading the antenna, has a dielectric constant of 30 instead of 58. The 
reason for this is that the supplier no longer has any materials with dielectric 
constants greater than 30. This has not created any serious problem since the 
exact value is not very critical. Voltage-standing-wave-ratlo (VSWR) measurements, 
performed on beef sanq>les and the calves of several human legs gave readings of 
1.5: 1 to 2.4: 1, depending on the thickness of the intervening fat layer. 

Since the radiometer is a broad band device, the antenna bandwidth was made 
slightly wider than that occiq>ied by the radiometer. This was accomplished by 
proper positioning of the end plate on the antenna and length of coupling probe. 

Due to the very low loss tangent (0.0005) of the dielectric loading material, 
the power loss through the antenna is almost negligible. 
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Waveguide Antenna 


Waveguide Brass 

Dielectric Ftycast K' • 30 

Emerson & Cummings Inc. 



3.4.2 BADIO FBfiQUENCY SECTION 


The RF section o£ the radiometer consists of two parts; 

1. The BF amplifier section lying between the Calibration Switch and 
Square -Law Detector, and 

2. The radiometer Calibration Section. 

All of the BF components have 50-ohm Input and output Impedances and are 
Interconnected with miniature seml-rlgld 50-ohm coaxial cable. In many cases, 
interconnections are accomplished directly with miniature SMA-type connectors 
satisfying the requlresients of MIL-C-39012. 

1. BF Amplifier Section 

Referring to Figure 3-3, after the signal from the antenna passes through 
the Calibration Switch, it is applied to the microwave diode Modulator Switch. 

This comptment is operated at a 1,000 Hz rate by a 50%-duty-cycle square wave 
generated by an oscillator in the Lock-In Amplifier. Thus, the switch dwells one- 
half the time on the signal and the other half of the time on a resistive termina- 
tion, coamonly referred to as the Reference Load. 

A Phase Adjustment circuit In the Lock-In Amplifier permits correct phasing 
of the 1,000 Hz square wave at the Modulator Switch so that its operation is 
synchronized with that of the coherent detector in the Lock-In Amplifier. 

The switch has an Insertion loss of 0.7 db, a VSWR of 1.55 : 1 and isolation, 
between terminals, of 70 db. 

An isolator is located between the Modulator Switch and the low-noise RF 
Pre-Aiq>lifier to provide isolation, for the pre-asplifier, from ispeda*' .« fluctua- 
tions at the output of the Modulator Switch. The isolator provides 23 db isolation 
and has a forward Insertion loss of 0.5 db. 

The BF portion of the radiometer is, basically, a sensitive fixed-tuned radio 
frequency (TBP) receiver, it was not necessary to resort to a superheterodyne con- 
figuration, in this design, due to the availability of low-noise solid stats RF 
amplifiers. 
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BF signal asipliflcatlon and bandpass shaping is provided by two low^nolse 
transistor aaplifiers and a bandpass filter. The first amplifier has a noise 
figure of 2.4 db and gain of 23 db. In the absence of the bandpass filter, the 
simplifiers exhibit a rather broad 3-db bandwidth, from 590 to 830 MHz (240 MHz). 
Although a broad operating bandwidth is desirable in a radiometer, fcr reasons 
of ii^roved tei^rature sensitivity, the lower end of this bandpass overlaps the 
IIHP television broadcasting band. Thus, to eliminate TV interference, a steep> 
sided bandpass filter was added to the amplifier section. The filter 3-db band- 
pass extends from 755 to 904 MHz (150 MHz); however, since the upper frequency 
cutoff, in the BF amplifiers, is at 850 MHz the net passband is 75 MHz (755 - 830 
Mis). The bandpass filter provides 61 db rejection at 716 MHz (Channel 54) which 
is the highest UHF TV frequency in the MASA-Anes area. Tests at NASA-Ames have 
shown this amount of rejection to be adequate for proper operation of the radio- 
meter. The insertion loss of the filter is 0.8 db, from 760 to 901 MHz. 

It should be mentioned that, at t*'e Badiometric Technology, Inc., facility 
a different bandpass filter and a band rejection filter were required to eliminate 
UHF TV interference. The bandpass filter allowed a 3-db bandwidth of 119 MHz, in 
the radiometer after rejection of Channel 56 with an Interdigital-type band rejec- 
tion filter. The characteristics of all three filters are given in Table 3-1. 

Table 3-1 

BADIOMETER FILTER CHARACTERISTICS 


Filter 

Type 

3-db 

Bandpass 

(MHz) 

Stop Band 
Rejection 

Insertion 

Loss 

VSUR 

Bandpass 


61 db e 716 MHz* 

0.8 db from 



(At NASA-Ames) 

755 - 904 

60 db e 936 MHz 

760 to 901 MHz 

1.5 r 

1 

Bandpass (At Radi 

.0- 

52 db a 660 MHz 

0.8 db from 



metric Tech., Inc 

:.) 700-890 

>70 db @ 625 MHz 

720 to 805 MHz 

1.5 r 

1 

Band Rejection 


3 db from 


1.5 r 

1 

(At Radiometric 


720 to 730 MHz. 




Technology, Inc.) 


60 db from 






723 to 728 MHz. 





^Higbaat interfering UHF TV frequency. 
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It ahould be nentloned that, although WF TV interference has been eliminated 
through the use of the above filters, some Interference is caused by fluorescent 
ll^ts. This is only present when the radiometer is being used in the normal 
Operate mode and amounts to approxim^ely 12°K at the output of the radiometer. 

Since this interference is relatively broad band and cannot be eliminated 
by filtering, the radiometer should be operated only in areas where incandescent 
lighting is available. Wherever possible, nearby fluorescent lights, in adjoining 
rooms, should be switched off. 

The second BF ao^lifler has a gain of 43 db and a noise figure of 4.1 db in 
the passband used at NASA-Anes. The total RF gain, represented by the two RF 
amplifiers and the intervening bandpass filter, is 65 db. The signal level at the 
input to the Square-Law Detector/Video Amplifier unit may be determined by calcula- 
ting the noise temperature of the receiver and, then, substituting the result in 
the equation, 

Pj, • kTB (3-2) 

where, 

PfI • equivalent noise power 
k ■ Boltsmann*s constant 
and B • receiver bandwidth 

The noise temperature of the RF amplifier section, referred to the input 
terminals of the first amplifier, may be calculated with the aid of the expression, 

NT ■ (F-l)T (3-3) 

where, 

NT is the noise cenqierature , °K 
F is the amplifier noise figure 
and T is the ambient tesiperature of the amplifier, °K. 
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Now F ■ 2.4 db (1.74 numerical value) 
and T • 3ll°K 

Therefore, NT ■ (1.74-1) 511 • 230®K 

Now, k ■ 1.38 X 10’^^ 

and B ■ 75 X 10® 


Thus, the equivalent noise power at the input to the first RF amplifier, is 

Pj, • 1.38 X 10"*’ X 230 X 75 X 10® 

■ 2,38 X 10"^’ watts 

or Pjj • 2,380 X 10"^’ milliwatts 


Converting to dbm, where db ■ 10 log 

P„ ■ 10(-13) ♦ 10 log 2,380 

II 



• -150 ♦ (10 X 3.376) 


or Pj^ ■ - 96 dbm. 


With a total gain of 65 db in the amplifier section, the signal level at the input 
to the Square-Law Detector is -31 dbm. 


luare-Law Detector/Video Amplifier 


The function of the Square-Law Detector/Video Amplifier unit is to 
recover the low frequency cosiponents of the rcdioaeter input signal and provide 
output voltages that are proportional to input power. The tangential signal 
sensitivity of the unit is -51 dbm. Thus, the -31 dbm signal level, calculated in 
the preceding paragraph, is adequate for proper operation of the detector. This 
input level does not overload the detector/asqplifier unit. Other characteristics 
of this unit are as follows; 


BF frequency band; 600 - 800 MBs 
Video Bandwidth; > 100 KHs 
VSWR; 1.7 ; I 
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3.4.3 C4UBR4T10N SECTION 


The purpose of this portion of the radiometer is to provide a means of 
establishing a zero baseline for all signal levels and to furnish a noise tempera- 
ture calibration signal, of known amplitude, so that output voltages may be 
accurately converted to noise temperature, at the input to the radiometer, at any 
point in a given time interval. 

As sho«m in Figure 3-3, the Calibration Section consists, in part, of a 
resistive termination (Base Load) connected to the input of the Modulator Switch 
via a directional coupler and the Operate/Calibrate coaxial switch. When the 
switch is in the Calibrate position, and the L'oise Source is switched off, the 
Modulator switch is operating between two resistive terminations at essentially 
the sane temperature - the Base Load and Reference Load - and the radiometer input 
is at a null. This establishes the sero baseline for all signal and calibration 
levels at the input of the radiometer. 

When the noise diode is switched on, it generates a noise temperature, in 
excess of the standard ambient temperature (290 degrees Kelvin) by a specified 
number of decibels. Thus, 

Excess noise temp. ■ db excess noise (numerical) X 290, (3-4) 

In this particular radiometer, the excess noise capability of the noise diode 
is 26.65 db, at a frequency of 800 MHz. This corresponds to a numerical value of 
462.4. Thus, in this case. 

Excess noise tenq>. • 462.4 X 290 ■ 134,096% 

To determine the noise tesqperature at the input to the Modulator switch, it 
is necessary to take into account signal attenuation, by the intervening lossy com- 
ponents, and re-radiation by the components at their ambient tesperature. This is 
accoiq>lished through the use of the well-known expression, 

L “ 
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or. 


■ T ♦ 


(3-6) 


No 



o 


K 


where , 

Tjj^ is the noise temperature at the output of the lossy componeiit 
Is the noise temperature at the input of the lossy component 
T is the thermodynamic temperature of the component 
and L is the attenuating loss through the component. 

The lossy components consist of the fixed attenuator, interconnecting coaxial 
cables and the Operate/Calibrate coaxial switch. The directional coiq>ler involves 
a 20 ± 1 db coupling loss, with negligible insertion loss. Listed below are the 
attenuation values of the lossy components; 

Fixed attenuator 16 ^ 0.5 db 

Interconnections and 

coaxial switch 0.5 db (estimate) 

Using Equation (3-5) iterative calculations \tere performed to detenaine the 
output noise temperatues of the above components, allowing a 20-db nominal coupling 
loss between them. The result is a noise temperature of 55^K at the input to the 
Modulator Switch. As discussed in Section 3.5, the value derived from over-all 
radiometer calibration, using a standard hot-cold noise generator, was 53.1°K. The 
agreement is rather close between the calculated and derived values. This indicates 
that the error Involved in using nominal losses for the attenuator and directional 
co<q>ler, and an estimate for the remaining losses, was fairly small. 

The above method of generating an internal calibration signal is simple and 
reliable, particularly if the ambient ten^erature of the coaq>onents is held within 
close limits over the periods of time involved in radiometric measurements. As 
stated in Section 3.3.2, the Br head temperature was stabilised at 311 i 0.5*^ during 
measurements performed with the radiometer. 

It is worth observing that good long-term stability and absolute accuracy 
are dependent on close control of RF head temperature. The above uncertainly of 
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t 0,5^K can be largely reduced through the use of a proportional temperature con- 
troller with finer resolution and the application of a thicker layer of insulation 
on the inner walls of the RF head. Calibrati^.. signal stabilitv can a'so be 
improved by keeping the DC power on the noise diode continuously and placing an 
on-off switch between the fixed attenuator and directional coupler. 

3.4.4 LOCK-IN AMPUt lER 

The Lock-In Amplifier is a standard fixed frequency commercial unit furnished 
by Princeton Applied Research Corporation. Their Model 120 is used in the radio- 
metric system. 

Lock-in amplifiers are used widely in radiometry due to their ability to 
measure very small changes in amplitude in slowly varying AC signals, which are 
usually obscured by high amplitude noise. In Equation (3-1) the noise temperature 
of a radiometer is expressed by the term (FL-I)Tq. The value of this term is 588°R 
In this radiometer. Wito 10-second output time constant, the radiometer is capable 
of measuring changes of 0.06*%. This represents 0.01% of the noise temperature of 
the radiometer. Clearly, the lock-in amplifier performs a very important function 
in the over-all system. 

This unit provides ten (10) selectable output time constants, from 1 msec, 
to 30 seconds and a connection for external low pass filters. It also furnishes 
the 1,000 Hz 50%-duty-cycle square wave, and a phasing control, for operation of the 
radioBieter input signal Modulator. The ± 45 degrees of phase shift, available in 
the phasing control, permits easy synchronization of the Modulator Switch and 
coherent detector located in the Lock-In Amplifier. 

The full scale sensitivity of the Lock-In Amplifier is 100 microvolts. It 
has a 24-bour stability of ^ 0.2% of full scale. 

3.4.5 RADIOMETRIC DATA SYSTEM 

The Radiometric Data Acquisition System (RDS) was designed tc acquire, 
digitize and record data from the microwave radiometer. Data is stored on conven- 
tional Phillips magnetic tape cassettes. 
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The electronics and the cassette recorder ar« housed i.a a portable cabinet. 

The eyetem can output Its data to a local teletype or to a remote computer using 
a Bell 202C data set. 

This Section describes the system, Its accuracy, data storage cai>ac'ty at.d 
modes of operation. 

The BDS has been designed with modern Integrated circuits using 7400 secies 
TTL with some medium scale Integrated circuits to minimize the component count. 

All units are standard plug-ln components to facilitate field repair when necessary. 

All data Is recorded In a serial digital format which confoimu to the ASCII 
standard for asynchronous data. 

Each character Is recorded serially on an lncrementt.1 cassette recorder which 
always records the bits of each character at its maximum bit rate. The model 133 
tape recorder accepts a 300-foot certified Phillips cassette and records both the 
data and the data complement on two track* fr>r . lundancy. 

The BDS system block diagram Is shown In Figure 3-5, which shows the inter- 
connections of all major components. This description will explain the Interaction 
of the major components. 

Badlometer - The radiometer outputs a signal which has a full scale range of zero 
to ^10 volts but uses a smaller portion of that range during actual data 
acquisition. An output variation of 0.1^ corresponds to a 10 millivolt (mv) 
change In output voltage. 

The BDS is designed to have an absolute accuracy of t 2 mv and a resolution 
of 1 mv. A peak-to-peak noise output of 60 mv must be averaged to detect the 
desired O.l^K changes In radiometer output. 

Analog to Digital Converter f\DC) - A 16-blt decimal, Integrating ADC was selected. 
The unit employs very stable voltagc-to-pulse count conversion with a basic 
accuracy of 0.01%. In order to Integrate Interfering 60 Hz noise, the unit 
Integrates each data sample for 16.7 ms. 
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ThumbwhMl 

8«l«etion 






Th« ADC has a hi^ iapedance dlCCarcntial aapllfiar which is preceded 
bp an operatlMial anpliCier, used for voltage scaling to accoemodate the 0 
to ♦lO volt radioaeter output signal. Nonully, the radioaeter is connected 
to the ADC and is constantly being digitized at one of the four rates (1, 2. 

4 or 5 sasples per second) selected by a rotary switch located on :!e oain- 
tenance panel. An operator can connect the radioaeter, energise the BD6 and 
aake any adjustaents that are necessary while wati^ing the ADC output. 

Analog Multiplexer - Solid state differential aultiplexing is used in the RD8 to 
connect thermistors and calibration voltages to the ADC. A 16-diannel analog 
anltiplexer is used with a 2-channel aultiplexer; the latter is shoun as a 
circular switching symbol in the block diagram. Whenever the RD8 user selects 
Tenperature (T) (ref. T^T line to system control) port #1 is connected to 
port #2 and the Analog MUX (16-channel unit) is connected to the ADC. 

System control advances the Analog MIK throu^ its 16 positions in 
response to a start consuoid when (T) is selected, or one step at a tine in the 
TEST node. The section labelled Analog MUX also contains the constant current 
generator that forces 257.9 aicroasperes through the thermistors ns each one 
is connected. Each thermistor is nominally 1,471 (dins at 35*^C, which is 
displayed as 0.350 volts on the 4-digit display. 

MIX Select logic - The Mux Select Logic (MSI) cerates under the direction of Systea 
Control (SC) to deliver the correct binary address to both the analog and 
digital aultiplexers. Logic within the multiplexers acconpliah decoding from 
binary to 8 lines or 16 lines. 

System Control - The system control (SC) logic contains the logic microprogram to 
operate the systM in the following aunner. 

In response to a start conmand, the EDS records the eight digital words 
inputted to the digital multiplexer as 16 ASCII characters. If T is selected 
the 16 analog MUX inputs are also digitised and recorded.* SC then proceeds 


*teader and T data are recorded at 60 diaracters per second regardless of the 
sampling rata selected. 
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to rocord radioueter sanples at the selected saapllng rate until the nuBber 
CO-99) of saaples specified by the (R) thunbiiheele have been digitised. SC 
finally ends the data sequence and lights the SDY light. 


If slope (S) is selected, the RDY light is not lit as long as the radio* 
■eter input is varying by note than the slope linit (ref. SET SIAPE). As soon 
as the radioneter has stabilised (NOSIAPE) SC limits BDY and the sequence can 
be started. The SC uses a precision 600 KHs crystal oscillator which is con- 
tained within the ADC unit as its basic logic timing reference. 

SC accepts digital saqiling rate CD6R) and output data rate (ODR) signals 
from the timing section. It uses these signals to derive the correct diaracter 
writing rate for the character formatter and the conversion rate which is one- 
fourth of the former (i.e. there are four characters generated for each sa^>le 
conversion). 

SC receives an inhibit signal from the maintenance pannl when the recorder 
is being tested (BECTEST). This causes SC to inhibit its normal cycle when 
STABT is depressed. SC creates a signal called "replace caosette** tdiich it 
derives from the situation where the HD8 ia in the record SNide and End of Tape 
CEOY) is sensed. The operator recognises this condition by obaerving EOT with 
the tape recorder in the record mode and takes the appropriate action by 
replacing the cassette. 

Write CT^aracter Formatter - The write character formatter QfCF) accepts 4-bit BCD 
diaracters and adds the correct 5, 6 and 7th bits to make eadi diaracter con- 
form to ASCII format. Each diaracter received by tte recorder is then 
converted to an 8-bit word with odd parity by the R133 recorder. 

Whenever thm ADC is presented with an overscale input (>9.99v) it records 
a question nark for the moat significant digit (TXXX). When receiving data 
from the digital multiplexer it accepts 8-bit data words and demultiplexes 
them into four bit bytes, adds the correct three most significant bits and 
records the resulting 7-bit word. 
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CMMttc Bftcorder - An increnental, serial cassette recorder Model #133 Banufactured 
by MsBodyne is used as the data storage slenent. Its salient characteristics 


arer 

Tape length 
Bit Packing Density 
Storage Capacity 
ReadAlrite speed 
SIBRT/STOP Tiae 
No. Tracks 
Record Fonuit 
Error Rate 


300 feet 
240 Bpi 

72,000 - 8>bit characters 
60 d>ar/aec. 

3 SIS and 7 as (reap) 

TWo 

NRZI - redundant 
1 bit in 10^ 


All conaands to, and data Coraatting for, the R1S3 are accosiplishsd in the RDS 
logic. The only function that an operator can perfora directly on the transport 
is rewind and this is done autoaatically when the cassette is inserted . 

Csssetta Becorder Control (GBC) the CBC is the control interface to the R133 and 
operates with the character fonuittcr to control the recorder and to output 
data. ORC conditioM the E/BOT and RW signals froa the tape recorder and the 
clean ti^ (CT) signal derived froa the read character fomatter. 


Bead Character Fomatter (RCF) • The &CF accepts parallel characters froa the R155 
recorder which asseables each group of B bits into a parallel -output. RCF then 
adds a start and a stop bit to the character and outputs the data as a serial, 
BS-252 coapatible, diaracter to either a Bell 202C aodea or a teletype at one 
of the selected playbadc rates. 


Although the playback character rates are 10, 30 and 60, their correspond- 
ing baud rates are not unifom aultiplea because the teletype (aodel ASR-33) 
requires a 110 baud output. The rate requireasnts are listed below. 


Character Bate 

Rate 1 10 Hs 
Rate 2 30 Bs 
Bate S 60 Bs 


Baud Rate 
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110 

300 

600 



Bate* 2 and 3 involve 10-bit characters with one start and one atop bit while 
Bate 1 requires 2 stiq> bits. While the 300 and 600 character per second rates 
are derived from a crystal oscillator, the 110 baud rate is produced by a 
snail auxiliary oscillator as shown in the block diagran. 

CleMi Tape lagic - When the SDS is used to playback data from a cassette that is 
not completely full, it is desirable to terminate the transmission as soon as 
possible to minimize playback tine. In order to do this the RCF has "clean 
tape** logic which operates in the playback mode and detects the fact that the 
playback head is now looking at clean tape. 

The clean tape logic uses the fact that a read start pulse is not followed 
by tape clock signals when the ti^ is clean. At any playback speed the tape 
clock poises will occur within 3 milliseconds of the start pulse so tbs logic 
for clean tape isr 

C1£AH XAFC ■ TAPE ChOCR for Apprtot. 10 ns 

The effect of clean tape (CT> is to extinguish the run mode by clearing the 
run flip-fl( 9 . 

Msintenancs Panel (MP> - The MP contains the necessary controls and displays to by- 
pass the BD8 front-end in order to record a test character on tape. There are 
also eight data lights to display characters that can be read one at a time 
while moving forward or backward. 

The recor^r test switch SEC TEST contained on the HP is used to place 
the system in that mode. This control and the sampling rate selector are 
really system function switdiaa but are located on the mabitananca panel to 
simplify the EDS front panal. 

Timing Section <1S) - The TS contains the necessary frequency divider logic to reduce 
the 600 Us reference clock signals to pulse rates which are i^plicable to the 
BD6 c^ration. Sampling rate selections are made by salaeting rates which 
exist in the TS logic. The TS logic also transmits the correct data sampling 
rate (D8R) and r^tput data rates (QOR) to the SC. 
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Data Storaf* Capacity - The BDS can record 72,000 aaynchronoua charactara on one 
300-foot caeaette. Since a fixed-length header ia recorded prior to each 
data record, and the nunher of temperaturea and calibration pointa are fixed 
when aelected, an equation can be develt^d to deacribe the BDS data atorage 
capacity. 


A convenient parameter la the number of data blocka where one block la 
defined aa all of the data between two headera. 


No. of Data Blocka 


■ 


72 X IC^ charactera 
N 

16B ♦ 4B^ ♦ 64T„ 

na 1 


In the above expreaaion, the variable B^^ can be 1 to 99 for eedi value of n. 

The value of R ia determined by two thumbwheel awitchea labeled B that are 
mounted on the BOS front panel. A aingle toggle awitch Mrked T and T deter- 
mlnea the value of Tn aa 0 or 1. The minimum number of data blocka that can 
be recorded ia given by the caae where the uaer alwaye eelecta tempkind 99 
redicawter aa^>lea. A *%iore- likely** caae ia where the uaer eelecta the tempera- 
ture (T) data ten percent of the time «»d recorda 20 data pointa. Both valuea 
are calculated aa followa; 


(Ho. of Blocka)^ 


72 X 10^ ^ 72 X 10^ 

■<ld ♦ 996 ♦ «4) 476H 


The largeat integer value for H • 151 which givea; 


(No. of Data BlockSj^i^^) • 
(No, of Data Blocka) • 

•Vw 

(No. of Data Blocka « 

ave 


151 

72 X 10^ 

16N « (64) ♦ 4 (20) » 

10 

72 X 10^ 

102 .4N 

700 
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5.5 BADIOMBTBR CAUBRA.TION 

The purpose of the radiometer calibration was to determine the relationship 
between output voltage and Input noise temperature, to assess the linearity of the 
radiometer response, to calibrate the internal noise tube, and to measure radiometer 
sensitivity. A hot>cold standard noise generator was used in the calibration. The 
noise ^nerator uses two resistive elements, one immersed in liquid nitrogen (77.3%) 
and the other mounted in a temperature-controlled oven (373.2%). The output of the 
noise generator was connected to a variable attenuator via an RF cable. The antenna 
was removed and an RF cable from the variable attenuator was connected directly Into 
the radiometer. The radlooieter input brightness temperature, Tg, can be cooq>uted 
using the radiative transfer equation, 

Where, 

T^^ * resistive load input brightness temperature (%) 
a resistive load temperature (%) 

T^^^ » resistive load output brightness temperature (%) 

Y ■ resistive load attenuation (db) 

Initially, setting T^q equal to the brightness temperature of the noise generator and, 
given the temperatures and attenuations of the three loads, consisting of two RF 
cables and variable attenuator. Equation (3-7) can be solved iteratively for the three 
loads to yield Tg. 

The calibration procedure was to record, 

1) The temperatures of the two RF cables and variable attenuator 

2) The variable attenuator setting and 

3) A series of 100 radiometer output readings (measured in volts) at 
a rate of one reading per second. 
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An avaragt and raa variation oC the 100 readings was then found. Peri* 
odlcally, the radiometer was switched to the Zero and Calibrate modes. In the 
Zero mode the radiometer was looking at the internal base load, and in the calibrate 
mode it was looking at the internal noise tube. As in the operate mode, 100 read- 
inga mere taken and V*, V^, and were found. Figure 5-6 shows (V^-Vp) 

plotted against Tg which was computed from Equation (3-7). Bach point corresponds 
to a different variable attenuator setting for which was measured. was 
essentially constant, varying by less than 1%. The radiometer response is very 
linear. Since we were particularly interested in the radiometer response at the 
higher temperatures, a linear fit, based on a least-squares appromimation, was found 
using the nind points corresponding to the nine highest Tg*s. The fit was, 

Vo-V^ ■ 0.05572 - 16.57 ± 0.01 mv (5-8) 

where, t o.Ol mv is the rms deviation from a straight line. Noting that the range 
of Tg over the nine points was 70 or, the uncertainty in the slope of the line is 
0.01 mv/70°K. The gain, G, of the radiometer is given by the reciprocal of the slope: 

G • 18.62 t 0.05 OK/mv (3-9) 

thm difference (Vc-V^) during the calibration was 2.65 mv. Multiplying this figure 
by the gain gives the amplitude of the pulse from the noise tube in <>K. The noise 
tubs calibration pulse was 53.1 ± 0.1%. Assuming this pulse stays constant, the 
figure can be used in subsequent experiments to determine Tg by means of the equation. 


where, 

is the temperature of the base load. 

The sensitivity, ATg, of the radiometer is defined as the rms fluctuation of 
the radicnster output, in the absence of an input signal, in % rms. 
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AX 0 ■ 6 AVo, ®K rma 


(5-U) 


ATgVas found to be close to 0.05^ in the calibration experiment, in which a 10 - 
■econd time constant was used. 

5.6 DATA SEDUCTION 

The radiometric data acquisition system (BDS) acquires and digitizes data on 
magnetic tape using a conventional Phillips data cassette . All data is recorded 
in a serial digital format which confonss to the ASCII standard for asynchronous 
data. The data is grouped into 'bund*» with carriage return and line feed characters 
separating the runs. Each run consists of the following data; 

CD Number of radiometer samples (l-lOO) 

(2) Experiment number (00-99) 

(3) Sampling rate (1, 2, 4, 5 samplea/aec) 

(4) Status (specifies radiometer mode and if temperatures were recorded) 

(5) A 4-digit number (for specifying time, attenuation, or other perti- 
nent information) 

( 6 ) Thirteen temperature probe measurements (if recorded) 

(7) Series of radiometer samples. 

The processing of this data was performed through National CSS, a time-sharing 
computer network using five m System 360 Model 67 processors. 

The existing MCSS software PTPtI£, used for Inputing asyndironoua data, 
required that a carriage return and line feed occur after every 160 characters, or 
laaa. This constraint required that the data recorded on the cassette first be 
editted. Because the cassette could not be stopped in playback node, a copy of the 
data was generated on punched paper tape by playing the cassette into a model 33 
taletypewriter. This punched pup^r tape was then played back into the teletypewriter 
to generate a sacond punched paper tape. During this second playback, the paper 
tape reader was interrupted after every 100 characters or so, and carriage return 
and line feed were manually typed in. These two additional characters appeared on 
the second punched paper tape, making it coi^tible with PTPllB. 
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Th« editte' paper tape was read into the computer on file IMPUTl. INPUTl 
eonaiated of a aeries of records, each record corresponding to the data that was 
sandwiched between sets of carriage return and line feed characters. Using an 
MOBS edit comaand, the character *'!'* was inserted at the end of the date in each 
record. The Fortran program SETUPl was then run. SETUPl reads IMPUTl, checks Cor 
errors in the formatting of the data, and blocks the data into separate runs on 
file IMPUT2. File INPUTS is then read by Fortran program SETUPS «diich checks for 
bad diaractera in the data field. Any bad characters are deleted and the remaining 
date is read onto file INPUTS. 

Fortran program RADIGM, along with associated subprograms BADBMS and PBTEMP, 
then read INPUTS. Thirteen titles, each consisting of 40 characters uk less, are 
neat inputed from the teletypewriter console into RADIGM. These titles serve to 
Identify the thirteen temperature probes. RAOIOM converts the temperature probe 
voltages to **C, finds the mean and rms values of the series of radiometer samples, 
determines how well the sanples agree with a random process, using the chl-sguare 
test, and uses Equations (3-10) and (S-11) to compute the apparent tei^rature 
and rms noise AT^, based on previous sero and calibrate ouide readily. This infor- 
mation, along with items 1 through 5 listed on the previous page, is printed out as 
shown in Figure S-7 . RADIOM also outputs onto file OS the probe tei^ratures 
and, onto file 04, the mean and rma values of the radiometer saaples and the chi- 
aquare test result. Listings of programs SETUPl, SETUP2, RADIOM. and sut^rogram 
PBIEMP and RADBMS are in Appendix A. 
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• • • ROS READ OUT • • * 


experiment NUHRER « 7 

time a 04t?l 

RAOtOMEtKR IN OPERATE MODE 

temperatures recorded 

Sample rate ■ i/sfc 

NUMBER OF radiometer SAMPLES «100 


RDS calibration 

ZERO PROBE SCALE FULL SCALE 
2 23S2 9997 


25*6 C 
37*4 C 
36«4 C 
37.0 C 

24.3 C 
37.9 C 

37.9 C 

34.9 C 

37.3 C 

• 27.8 C 

• 2B.2 C 

• 27.5 C 

m 26.4 C 


MEAN RADIOMETER OUTPUT • 2942.8 
RMS DEVIATION « 14.1 

CN1«S0UARE TEST • 0.24JE 00 


APARENT TEMPERATURE « 305.S7 DEO K 

RMS NOISE *1 O.IS DFG K 

RASED ON calibration AT TIME « 04110 


PRORE temperatures (DEO C) 

R1 • ROOM AIR 

R2 • BASE LOAD 

R3 a rectal 

R4 m RIphTSURFACE 

RS a LEi I SURFACE 

R6 a right one INCH DEEP 

R7 a right ONE^HALF INCH DEEP 

RB a noise SOURCE INPUT PONT 

R9 ■ RIGHT SUBCUT. 

RIO m left one-half INCH DEEP 
Rll ■ left subcut, 

R12 • LEFT ONE INCH DEEP 
R13 a ROOM AIR 


Pigura 3>7 - Sample Computer Printout of Badiometric and T\tmparature Probe Data, 
for Fijial Dog Experiment 
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Section 4 


TEMPERATURE MEASUREMENTS ON DOGS 
4.1 MEASUREMENTS ON FIRST DOG 


An investigation was conducted on a large mongrel dog, at the New England 
Medical Center Hospital, to attesq[>t to determine the relationship between blood 
Clow and surface and subsurface temperatures. This work involved the radiometric 
data system (RDS) only, with a variety of temperature probes and a flowmeter; the 
radiometer proper was not employed in these initial observations. Dr. W. Carney, 
of the Hospital staff, assisted with this experiment. 

The dog was anesthetized and the following ten^erature probes were installed 
in the inner thi^ of each hind leg; surface, subcutaneous, 0.95 cm deep and 1.9 cm 
deep. The dog*s abdomen was opened and a Carolina Medical Electronics, Inc., square- 
wave electromagnetic flowmeter was attached to the right iliac artery, leading to 
the test leg. The right iliac was then clamped shu** and ten^ratures were recorded. 

The results are shown in Figures 4-1 and 4-2. It should be noted that the 
initial downward trend of temperatures in both legs ia mainly due to the cooling 
effect of the laboratory air conditioner; the room tesq>erature h«d risen to a level 
of 82°F. The important feature, in the four pair of plots, is the difference in 
temperature between the reference and test legs. Thus, after 0.5, 1.0 and 1.5 hours, 
the drops in the rest leg temperatures were as listed in Table 4-1. 

Table 4-1 

IBfFEBATURE DROPS IH 1SST ISG 


Probe Depths 

After 0.5 hr. 

After 1 hr. 

After 1.5 hrs. 

1.9 cm 

0.65°C 

0.70^ 

0.90OC 

0.95 cm 

0.60 

1.05 

1.50 

Subcut. 

1.70 

2.40 

3.20 

Surface 

0.80 

1.40 

2.00 
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The above data shows the largest drops in the subcutaneous probe, with the 
surface probe being the next largest. The drops at the 0.95 cm depth are not too 
far removed from those at the surface. It should be noted, however, that collateral 
circulation would tend to compensate, to a considerable extent, for blood cutoff in 
the iliac artery. Thus , the above teaperature drops are quite encouraging. 

It is difficult to explain the continued increase in test leg tenq>eratures, 
for the 0.95 cm and subcutaneous depths, after 11:45 am; the corresponding refer- 
ence leg temperatures were decreasing after this point in tine. It appears that 
the increased blood flow had a greater effect at these depths than at the surface 
and 1.9 cm depth. 

Unfortunately, the readings lumished by the electromagnetic flowmeter were 
somewhat erratic; thus, the blood flow in the right iliac could not be determined. 

4.2 MEASURBIENTS ON SECOND DOG 

A large mongrel dog, similar in size to a German Shepherd, was anesthetized 
and prepared for ni^ttime radiometric observations, after the IHF TV stations had 
gone off the air, to avoid TV interference with the radiometer. Temperature probes 
were installed in the thighs of the hind legs as follows: surface, subcutaneous, 

1.3 cm &nd 2.6 cm. Instead jf clamping off the right iliac, as in the experiment 
described in the above Paragraph, a tourniquet was applied to the upper thi^ of the 
right leg. It was expected that this would reduce collateral circulation and the 
flow in the arteries would also be smaller; it was also expected that a faster 
temperature ^aponse would be achieved. A Statham Model M-4001 Square-Wave Electro- 
magnetic Flowmeter was arranged for use in this experiment. Dr. W. Carney was in 
attendance during this investigation. 

Prior to applying the tourniquet, voltage standing wave ratios were measured 
with the antenna disconnected from the radiometer and placed on the dog's right leg. 
The average value was 2.6: 1, across the radioemter bandwidth. It can be shewn that, 
for this value, the antenna power transmission is 0.8 (see Section 5.2). 
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Following this procedure, the antenna was re>installed on the radioaeter and 
placed in close contact with a cleanly shaven area near the probes on the ri^ t leg. 
Teaperatures were, next, recorded on both legs and the tourniquet applied as L'.ghtly 
as possible to the upper thigh of the right leg. Radioaeter readings were taken and 
teaperatures recorded for a period of 50 ainutes. Figure 4-5 shows that the resultant 
radioaeter apparent teaperatures dropped 1.7^ during this period. Thermodynaaic 
teaperature changes, recorded during this period, are listed in Table 4-2. 


Table 4-2 

TEMPEBATOBE CB4MGES IH REFEREMCE AMD TEST LEGS 


Probe 

Ref. Leg 

Test Leg 

Depths 

Teap. Change 
(50 ain.) 

Teap. Change 
(50 ain.) 

2.6 ca 

- 0.5®C 

- 1.1®C 

1.5 ca 

- 0.4 

- 2.1 

Subcutaneous 

- 0.5 

- 1.7 

Surface 

Ho reading 
(probe fell off) 

- 2.1 

Rooa tMperature changer ♦ 0.4*^C 


Many aore readings were taken over a period of 5.5 hours; however, imfortu- 
nately, all of the data was lost due to nuaerous spurious characters recorded by 
the radioaetric data recording systM (RDS). This was caused by overheating of the 
aagnetic tape recorder. The above saaple data was recorded aanually froa the RDS 
display. Based on these results, it is probable that the r^iooeter sensed subsur- 
face te^>erature changes froa a depth of, at least, 1.5 ca. 

The flowaeter probe was applied to the feaoral artery and vein in the test 
leg, during the above observations. A rigid plastic tube was placed around these 
blood vessels and the tourniquet was placed over the tube and around the thigh of 
test leg. 
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Tourniquot Appliod Tourniquot Raloaaod 


Figure 4-4 eboHs the 'ouneter reading! versus tine. The tourniquet was 
applied at 5;25 am; however, the flownater did not show any drop In blood flow 
until the femoral artery was clamped shut at 6:25. When the clamp was opened S 
minutes later, the flowmeter readings gradually rose to a rather high value of 108 
ml^ln., and then began dropping to 66 ml/mln. , with no other changes made In the 
experiment setup. Due to the unpredictable nature of the raolngs, the flowmeter 
was not used In any furtter sieasurements . 

The discouraging results obtained with two different flowsmters suggest that 
other methods be explored for use In correlative blood flow measurements. It Is 
understood that radioisotope techniques could provide useful Information on blood 
flow for this type of experiment. 

4.5 ICASISEMERTS OM THISD DOG 

A large black Labrador male dog, weighing 57 lb., was used In this experiment, 
lbs dog was anesthetized with 11 cc of pentobarbital and four (4) temperature probes 
were attached to, and inserted In, the thighs of eadi hind leg. Dr. H.H. Miller was 
in attendance. The right hind leg was chosen to be the test leg. A temperature 
probe was Inserted In the esophagus. 

The radloamter was used in making measurements on the test leg; however, due 
to excessive UHF TV Interference, none of the radiometer data was of any value. 
Bo n et he leas, the extensive temperature readings, taken during the course of the 
experiment arm of some Interest. 

Figures 4-5 and 4-6 show temperature readings for various events throughout 
the 6-hour experiment. Referring to Figure 4-5, it will be noted that, during a 
56-minute interval, following clasy in g of the femoral artery In the teat leg, the 
only significant temperature drop occurred at the surface of the leg. After a tour- 
niquet was applied to the leg, there was a drop of 1.4^0 at depths of 1.3 cm and 

2.5 cm, over the next 59 minutes. It will be noted that the subcutaneous and sur- 
face te^teratures increased 0.6<\! during this period; this is probably because the 
room temperature rose by 2,5^C in this interval. 
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Following release of the tourniquet and femoral artery clamp at 5;23 am, the 
surface, 1.3 and 2.5 cm temperatures rose steadily by 3*^C over the next 32 minutes; 
however, the subcutaneous temperature rose only about 0.4^C during this time. It 
is difficult to explain the latter, except that the blood flow could have been 
relatively limited in the subcutaneous region, thus retarding a rise in temperature. 

Figure 4-6 shows a continuation of the temperature data, obtained under some- 
what different conditions, involving the use of a tourniquet, femoral artery clamp 
and icebag containing an ice-alcohol mixture. Unfortunately, this plot has a data 
gap in it, due to some equipment problems; however, the effect of the icebag is 
evidently quite pronounced at the surface and all depths between 7:23 and 7:47 am. 
The temperature reduction is greatest, of course, at the surface of the test leg 
and becomes progressively smaller at the subcutaneous, 1.3 cm and 2.5 cm depths, 
respectively. This is understandable since the effect of a cold body, placed on the 
surface, would tend to be reduced with increasing depth, due to the thermal insula- 
tion afforded by the intervening tissue. 

When the icebag was removed, at 7:47 am, the surface and subcutaneous temp- 
eratures began to rise while the 1.3 and 2.5 cm temperatures continued to drop for 
a few minutes. This lag, in reversing the downward trend in temperature, is greater 
at the 2.5-cm depth. The temperature at this depth finally began to rise 19 minutes 
after the icebag was removed and 3 minutes after the tourniquet was removed. Thus, 
the restoration of normal temperatures at depths of 1.3 and 2.5 cm was relatively 
slow until full circulation was permitted by the removal of the tourniquet. This 
Indicates that very little, if any, collateral circulation was permitted by this 
device . 


The above data shows that an ice mass can be employed quite effectively in 
establishing conditions suitable for investigations of the heating effects of blood 
flow at different depths in the body. This technique can be readily and safely 
adapted to blood flow investigations in humans wherein a microwave radiometer is 
employed in sensing subsurface temperatures. 

Alternatively, a heat pulse could be applied, at the location of interest, 
with the radiometer measuring the time taken to cool that area. In this connection, 


4-11 



an interesting experiment is described in Reference 3. The technique involves 
applying a heat pulse to the skin and measuring the rate of blood How with the 
aid of two thermistor sensors placed a known distance apart. This method could 
furnish useful correlative data during radiometric observations. 
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Section 5 


RADIOMETRIC OBSERVATIONS 


5.1 INTRODUCTION 


As discussed in Section 4, early attempts at obtaining radiometric data, on 
phautom models and blood flow in dogs, were thwarted by excessive externally 
^nerated interference. To determine the exact frequencies of the interfering 
signals, a survey was conducted, in the laboratory and at the New England Medical 
Center Hospital, using a sensitive spectrum analyzer. Based on information obtained 
in this manner, suitable bandpass and band-re jection filters were specified and 
purchased for the radiometer. These included a bandpass filter for use at NASA- 
Ames. The characteristics of these components were described in Section 3. 

Following incorporation of the filters, the radiometer performed satisfacto- 
rily for phantom measurements and observations on dogs. 

5.2 PHANTOM MEASUREMENTS 

TWo phantom models were established and measurements were performed on them 
in the laboratory. Data obtained from the first model indicated an excessive amount 
of outside RF interference and no conclusive results were obtained. The second 
phantom, shown in Figure 5-1, was entirely enclosed in a double layer of aluminum 
foil to eliminate any possible interference. To prevent evaporation the beef fat 
and lean beef were also wrapped in Saran Wrap. Four temperature probes were 
inserted at the various depths shown. Radiometer observations were made concurrently 
with temperature probe measurements as current was gradiially applied to both heaters. 
A 10-second time constant was used, and 100 samples were taken, at a rate of 1 
sample/sec., for each radiometer observation. During a 30-minute period, shortly 
after the beginning of the experiment, a partial short in the heater resulted in 
the radiometer measurements being in error. 

The results of these measurements are shown in Figure 5-2. The highest 
temperature throughout most of the 2-hour period was shown by a probe at a depth of 
5.1 cm, which was the closest to the heaters. Enclosing the phantom in aluminum 
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Figure 5-1 - Phantom Model 
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Figure 5-2 Phantom Apparent Temperatures and Temperatures 
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foil resulted in moch o' the heat convecting up the sides to the top surface. This 
caused the surface temperature probe to be warmer, at times, than the 1.9 and 3.8 
cm-deep probes. The short-dashed curve in Figure 5-2 is the measured radiometer 
apparent temperature given by Equatior (3-10). ihe layer depths shown in Figure 
5-1 and the four probe temperatures, given in Figure 5-2, were inputted into pro- 
gram lAYER, discussed in Section 2, along with the appropriate dielectric permit- 
tivities; theoretical brightness temperatures, Tg, and reflectivities, R, were then 
computed. Equation (2-1) was then used to derive theoretical apparent temperatures, 
shown by the long-dashed curve in Figure 5-2. The agreement between the measured 
and theoretical apparent temperatures is excellent. Considering the uncertainties 
in specifying rmittivities and temperat«*re profiles in the layers, and the fact 
that the phantom did not have au infinite cross-section (in fact, the phantom cross- 
section is only a fraction of fre». space wavelength), the small 1®K offset 
between theory and measurement is remarkable. 

The effec«,ive temperatuie, Tg, of a layered medium can be defined by an 
equation analogous to Equation (2-1). 

T^ = (l-R)Te + RTp (5-1) 


or, 


Te = Tg/ (1 - R) (5-2) 

The reflectivity. R, can be determined by measuring the voltage standing wave rstio 
(VSWR) denoted by S. 


R 


1 


4S 


(S +1)^ 


C5-3) 


The transmission t is de^^ined as 

T = 1-R (5-4) 

and is, plotted versus S in Figure 5-3. 


5-4 



POWER TRANSMISSION 
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0.5 



VOLTAGE STANDING WAVE RATIO 

Figure 5-3 - Power Transmission versi<s Voltage Standing Wave Ratio 
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By aeMuring T^, and S, can be found. S was oeasured to be 2.0. 

At 6;0V pa, when the tenperature spread throughout the phantoa was the greatest, 
equalled 50.5*’C. is a weighted average of the tenperatures throughout the 
the layered aediias, with the vei^ting depending on the power penetration depth. 

The tenperature 50.5^C corresponds roughly to the tenperature at a depth of 3.8 cm. 
Thus, it can be stated with some confidence that the radiosmter sensed teuiMsioturcs 
to a depth of, at least, 3.8 cm and, possibly, down to a depth of 5.1 cm. 

5.3 tSASUKEMEMTS ON lARCX DOG 

Badiometer measurements were performed at the New England Nedical Center 
Hospital, on a Labrador male mongrel, weighing 57 lbs., with Dr. Harry H. Miller 
in attendance. The dog was anesthetized wita pentobarbital and the left rear thigh 
was closely shaved. Temperature probes were pos:tioned on the thigh surface, sub- 
cutaneous, and at depths of 1.3 and 2.5 cm. The dog and radiometer were placed in 
the X-ray room, to avoid interference generated by the fluorescent lights in the 
Surgical Research Laboratory. Other external interference was fairly slight. A 
10-second time constant was used in measuring apparent tem|.eratures, and 100 samples 
were taken at a rate of 1 sample/sec., for each measurement. 

The measured radiometer apparent temperature and the temperatures of the 
surface, subcutaneous, and 1.3 cm-d.?ep probes are plotted in Figure 5-4, Figure 5-5 
shows the measured radiometer apparent tesperature and the 2.5 cm deep probe plotted 
on an exaggerated vertical scale. The vertical bars, in Figure 3-5, indicate the 
measured radiometer sensitivity, AT^^. As indicated by the annotations on the 
figures, the procedure followed consisted of the following steps: 

1. Reference radiometer and temperature probe measurements were 
initially taken. 

2. The radiometer was removed at 2:00 pm and repositioned and more 
reference measurements were made. This allcwed for checking on how 
sensitive the measurements w -.re to repositioning. Tlie apparent 
temperature curve indicates that the repositioning did not have a 
significant effect on the measurements. 
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Figure 5«4 - Apparent Temperature and Temperature Data «i Dog's 
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Figure 5-5 - Apparent Temperature and Tei^rature Data on Dog's Thigh 
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3. The radiometer was removed at 2:30 pm and a tourniq iet and ice pack, 
containing an ice-alcohol mixture, was applied to the upper portion 
of the thigh. After the leg was cooled down, the ice pack was removed, 
the radiometer was repositioned, and Bieasureinents were repeated. As 
expected, the apparent temperature showed a sharp decrease. 

6. The radiometer was removed again at 3:05 pm, and the ice pack was 
reapplied for further cooling. After sufficient cooling, the ice 
pack was removed, and the radiometer was placed back on the thigh. 
Measurements tier's taken as the thigh warmed up, during which the 
apparent te^>erature showed a steady increase. 

5. At 4;05 pm the radiometer was removed so that the tourniquet could 
be released. The radiometer was then repositioned and the apparent 
temperatures continued to increase. 

6 . To aid in warming the leg, a hot water bottle was placed on the 
abdomen between 4:30 and 4;45 pm. 

7. At 5:00 pm the radiometer was removed, the hot water bottle was placed 
on the inside of the leg for several minutes and then removed, and the 
radiometer was repositioned. 

Referring to Figure 5-4, it will be noted that the apparent temperature curve 
follows the 1,3 cm-deep temperature curve fairly closely. In Figure 5-5, the 
correspondence between T^ and the 2.5 cm temperature curve is rather poor when the 
temperature is increasing except, perhaps, after 5:00 pm. Thus, on the increasing 
temperature portion, the apparent temperatures appear to follow temperatures at 
shallower depths, as indicated in Figure 5-4. It will be recalled that the penetra- 
tion depth in the phantom model was somewhat greater - between 3.8 and 5.1 <ai. This 
is probably due to the lack of circulating blood in the phantom. However, the above 
results are quite encouraging and it appears that indirect blood flow measurement, 
with a microwave radiometer, holds considerable promise. 
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Section 6 


CONCLUSIONS 


Several cmcluslons of importance may be drawn from the results of this study. 

These are as follows; 

1. The development of a highly sensitive microwave radiometer, for use 

in near-field measurements of radiation from biological materials, has 
been successfully accomplished. Such observations are new quite 
practical and can be made with a high level of absolute accuracy. 

2. Measurements made on a phantom model indicate that the radiometer is 
sensitive to temperature changes occurring In the top 0.5 cm layer of 
fat and, at least, the first 5.2 cm of lean beef. It Is likely that 
temperatures as deep as 5,1 cm Influenced the radiometer response. The 
ability of the radiometer to sense temperature changes, from these 
depths, establishes the feasibility of measuring subsurface temperature 
chants In biological tissue. 

3. The close agreestent between the phantom measurements and the correspond- 
ing theoretical apparent temperatures, ^substantiates the validity of 
using a layered tissue model, of Infinite lateral extent, for predicting 
brightness temperatures In real tissue structures. 

4. Measurements performed on a hind thigh of a dog Indicate that the 
radiometer sensed temperature changes from a depth of, at least, 1.3 cm. 
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Section 7 


RECOMMENDATIONS 

The results of this study indicate several areas in need of further investiga- 
tion, to properly develop the mlcro^iave radiometric blood flow sensing technique. 
These areas are presented below. 

7 A INTECStAL VSWR MEASUREMENT CIRCUIT 

A VSWR measureoient circuit should be designed and fabricated to mount inside 
the radiometer RF head. Tnis will permit convenient and rapid measurement of VSWR, 
via the radiometer antenna, during future observations on phantom tissue models, 
dogs and hiasan subjects. The use of bulky, standard laboratory equipment, for this 
purpose, is awkward and tio&e consuming. It also requires de-mounting the antenna 
from the RF head. The recommended integral circuit would eliminate these problems 
and provide VSWR information easily and accurately* 

7.2 MEASURBfENT BLOOD DISTRIBUHC^i IN THE BODY 


Due to problems caused by externally generated interference, the work in this 
study was sufficiently delayed to preclude measurements on human subjects. Accor - 
ingly, it is recomaended that an extensive investigation be performed on humans to 
develop a technique for determining blood distribution in the body* The radiometer 
could be employed in measuring apparent temperatures of the head and other parts of 
the body, under s treses imposed by heat and exercize. Of particular interest are 
heat losses from various parts of the body, under specified conditions* These are 
unknown at present. 

The radiometer is capable of sensing subsurface temperatures down to a depth 
of about 1.5 cm. Since a significant portion of the blood volume is in this region, 
measurements at different points or the body should yield a good indication of blood 
distribution and heat flow. 
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7.5 


CORRELATIVE MEASURQIENTS OF BLOCO FLOtf 


As discussed in Section 4, it has been difficult to obtain accurate cor- 
relative data on blood flow during temperature and radiometric measurements. Even 
if the square-wave electromagnetic flowmeters had performed satisfactorily, their 
use is limited to a single blood vessel; thus, no indication is given of the amount 
of collateral flow. The other disadvantage of these instruments is that a surgical 
operation is required to permit their use. This precludes their application in the 
type of investigation described in the above Paragraph. 

In view of the above problems, it is recoomended that correlative data on 
blood flow be obtained for future radiometric investigations by either a radioisotope 
method, an impedance plethysmograph, the ultrasonic method, a cold mass or a heat 
pulse (described in Section 4) whichever is most appropriate fcr a given experiment. 
If the necessary equipment, for a particular investigation, is not readily available, 
consideration should be given to performing the experiment at a facility that can 
furnish it temporarily, on a rental basis. 

7.4 LOWER FREQUENCY RADIOMETER 

The operating RF center frequency of the radiometer, delivered to NASA-Ames, 
is near 800 MHz. This is about 75 MHz higher than was originally intended and was 
caused by the need for a bandpass filter capable of eliminating certain UHF TV 
frequencies. 

Since the depth from which the radiometer can sense temperatures appears to 
be limited to about 1.5 cm, it would be worthwhile furnishing an additional RF head 
to operate at a lower frequency and, hence, provide a means for sensing temperatures 
at greater depths. Figure 2-8 shows that operation in the vicinity of 500 MHz should 
be optimum. The tendency of the antenna aperture to increase to objectionable 
proportions could be eliminated by raising the dielectric constant of the loading 
material to a value that would compensate for the increase in aperture. It has been 
recently discovered that loading materials with dielectric constants above 50 are 
manufactured by the Raytheon Company. 
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Operation at a lower frequency could permit temperature measurements on 
organs located at some depth in the body. An ideal approach would Involve a 
3-frequency radiometer. Since each frequency would be limited to a particular 
depth* such an instrument could be designed to three temperatures simultaneously* 
in profile* to a depth dictated by the lowest frequency. This approach would permit 
a better understanding of heat flow from the body. 
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APPENDIX A 

LISTINGS OF COMPUTEPv PROGRAMS 

SETUP 1 
SETUP 2 
RADlC»i 
PBIEMP 
RADRMS 



PAGF 1 


<;FTUP1 FOPTPAN pi Il) = »n| OH/U/73 IP.lA.lb 

VP/CSS ^4ATIO\AL C'iS. INC. (STamfOMP DATA CFNTER) 


niMFNSION A ( 130) •rMbfiO) «AThST (4) 

OATA ATFST/«< «*•> tt»= *f»? •/ 

data ATtSTX/*! •/ 

* data ATEST0/«0 •/ 

DATA ATtST?/»? •/ 

data ATEST3/»3 •/ 

OATA ATESTQ/'Q •/ 

PFAOdvlflOO) A 
1000 FORMAT (130A1 ) 

I COMTINUE 
IREG«1 
lENOsiaO 

IF(A(3) •FO,ATF.ST9,ANO,A(4) .F0.ATEST9) GO TO Q99 
10 CONTINUE 
JsO 

no ?0 t*lHEG«ItPO 

J«J«1 

1F(A(J) .FO.ATf SIX) GO TO .10 
IF(I.GT.SOO) GO TO 40 
ISAVE«I 
H(I)sA(J) 

20 continue 

30 CONTINUE 

RFA0( If 1000) A 

IF(A(11) •FQ.ATFST(4) .ANO.Ada) .rU.ATEST(4) .aNO.A( 13) .F0.ATCST(4) 
X.ANDf Ad4) .FOfATEST (4) ) GO TO SO 
1BF6*ISAVF*1 
ieN0«IRE6fl?9 
60 TO 10 
40 CONTINUE 

WRITE JStSOoo) 

5000 FORMAT (23m OATA IN INCORRECT FORM) 

CALL FXIT 
so CONTINUE 
1ST0R«ISAVE 
IRE6«1 
IEN0»14 

WRITE(2f200n) (H d I f I«lHEGf lEND) 

2000 FORMAT (14A1) 

IF(H(S) .EO.ATEST (i) .OR.H(S) .FQ.ATFSTI4)) GO TO 80 

on 60 I«lSf78 

TSAVE«I 

IF(8(T! .EO.ATFSTO.ANO.Hdfl) ,FQ.ATEST0,AN0.od4?) .EO.ATESTO 
X.ANO.flC I«4) .FO.ATFSTH.ANO.Hd^S) fFQ.ATEST3.AN0.H( I^fl) 
X.rO*ATFST9,ANO,HdfR) .FO«aTEST9.ANO.B(I«10) ,FQ.ATEST9I 
K60 TO TO 
60 CONTINUE 

WRITE (2f?0fll ) ATFSTOf ATFST9fATEST9f ATFST9 

2001 FORMAT (4A1) 

GO TO I 

70 CONTINUE 

IBEG-ISAVE-S2 

IFNO«ISAVE4ll 

WRITE(2*200?) (Bd) f IsIHEGf lENO) 

2002 FORMAT (64A1) 


SETQOOlO 

SETOOOPO 

SET00030 

SET00040 

SETOOOSO 

SET00060 

SET00070 

SETOOOSO 

SET00090 

SETOOlOO 

SETOOllO 

SET00120 

SET00130 

SET00140 

SETOOISO 

SET00160 

SET00170 

SETOOISO 

SETOOlOO 

SET00200 

SET00210 

SET00220 

SET00230 

SET00240 

SET002S0 

SET00260 

SET00270 

SET00280 

SET00290 

SET00300 

SET00310 

SET00320 

SET00330 

SET00340 

SET003S0 

SET00360 

SET00370 

SET003H0 

SET00390 

SET00400 

SET00410 

SET00420 

SET00430 

SET00440 

SET004S0 

SET00460 

SET00470 

SET00480 

SET00490 

SET00500 

SET00510 

SET0OS20 

SET00530 

SET00540 

SET00550 
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SFTiJPl FOHTRAn pi 


TI)*‘?TT 


OH/14/73 


1P.34.15 
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HO CONTINUE 
IBFfi«!£ND*l 
lENOsTSTOP 

NUMRAO* ( IFMn-lREG* 1 ) /4 
TPF6*If ND-4*NUMWAI»4l 
WRITE (P*P003) NUMfvAO 
?003 FORMAT (I3> 

WRITE (?*?00?) ,’H( I ) «lsIBEf>*!F.NU) 

60 TO 1 
909 continue 

WRITE rP* pool ) ATFST9#ATFST9f ATFST9*ATEST9 

call exit 

END 


SETOOjnO 

SET00570 

SET00560 

SET00590 

SET0060C 

SET00610 

SET006P0 

SET00630 

SET00640 

SET00650 

SET00660 

SETOC670 

SET00680 
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SFTUP? FORTHAN R1 TD=wTI OH/14/73 12.3h,08 PAGE I 

vp/r«;s nati<»nal css* inc, (stamforh data center) 



DIMENSION P(SOn)«A(lfl) 

gETOOOIO 


data AX/»* •/ 

SET00020 


DATA AY/»? •/ 

SET0003^ 


DATA A/M ♦*»? •••) »««4 *,*S «, 

''.ETOOOAO 


X»6 •*•? »«*« •*!<) •«•() •/ 

SETOOOSO 

1 

CONTINUE 

SET00060 


RFAO(1*1000) (M(I)«Isl«14) 

SET00070 


WPITE<?*1000) (HU) tlsltU) 

SET00080 

1000 

FORMAT (UAl) 

SET00090 


IF(H(3) .FO.ACN) .ANO.H(4) •K).A(N) ) CALI EXIT 

SET TOlOO 


IF(B(S) .EQ.AX,OR,ri(S) .FQ.AY) G() TO 35 

SETOOllO 


PFADUtlOOl) (H( I) «Isl«h4) 

SET0C120 

1001 

FORMAT (64A1 ) 

SET00130 


IF(R(li.EO,A(9) .ANO.Hi?) .FQ«A(9) .ANn.H(3) .FO^AiR) 

SET00140 


X.AND*H(A) •FO.A(R) ) GO TO RR9 

SET00150 


DO 30 T«1*6A 

SET00160 


DO 10 JsltlO 

SET00170 


IF(R(I) .EQ.A(J) ) GO To 30 

SET00180 

10 

CONTINUE 

SET00190 


KRF6«( (I-l)/4)*4*l 

SET00200 


KEN0»KRE643 

SHT00210 


DO ?0 K*KHFG*KFNU 

SET00220 

?0 

B(K)sA(S) 

SET00230 

30 

CONTINUE 

SET00240 


WPITE(2*1001 ) (H(l)*Isl«E ) 

SET00250 

3S 

CONTINUE 

SCTU0260 


RFADUtlOO?) NUMHaD 

SET00270 

100? 

FORMAT (13) 

SET002B0 


IEN0«NUMRA0«4 

SET00290 


RFAD(1«1001) (M(I) tIsl.iEND) 

SET00300 


IF(NUMRAD,LT,3) go T' «0 

SET00310 


NCOUNT«0 

SETti0320 


00 60 I«5*IFNO 

SET00330 


DO SO JBltlO 

SET00340 


IF(R(I) •FQ.A(U) ) GO To 60 

SETS0350 

**0 

continue 

SET00360 


ncount«ncount«i 

SET00370 


IF(NC0UNT,F0*2) Go ”0 RRfl 

SET 00 330 


K*4 

SFT00390 


IF(IFN0-I.l T*4) K*-4 

SET004C0 


R(I)BS(I»K) 

SET00410 

60 

CONTINUE 

SET00420 


SO TO 100 

SET00430 

SO 

S(S)«A(10) 

SET00440 


B(6)«A(10) 

SET004SO 


B(7)bA(10) 

SET00460 


R(S)«A(10) 

SET00470 

100 

continue 

SET00480 


BIDbRIS) 

SET00490 


B(?)bH(6) 

SETOOSOO 


B(3)bBC. ' 

SETCOSIO 


fl(4)BB(H) 

SE7(;G320 


WR1TE(2*1002) NUMMAO 

SET00S30 


WRTTE(2*1001) (B(I)*Ist*lENO) 

SET00S40 


60 TO 1 

SET00550 
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08/14/73 


12.38.08 
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9«8 

NU*H)Ans994 


SETOOS60 


WRITE(2»100?) 

NUMmaO 

SET00S70 


GO TO 1 


SET005B0 

9 <)<) 

PRITE(2tl00U 

(8m «Ultft4l 

SET00590 


GO TO 1 


SET00600 


END 


SET00610 
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PAOTOM FOPTPAM PI IIJSPTI OH/U/73 1?*3T«H PAGf 1 

VP/CSS — - NATIONAI. CS«;t INC. (STAmFOPD hat a CFNTER) 


DIMENSION TITlM< 10»13) 
rilMENSION ATFST<») 

DIMENSION IVOL T(100)«TFmPP(13) 

DIMENSION CALRC?) vITImFCA) 

data ATEST/»< •••> •••* •••? •••> »t 

n*% 't*? */ 

ITIMEC1)»0 

ITIME|?)«0 

ITTME|3»«0 

ITIMEC4>«0 

CALRU^sO* 

CALR(?)*U 

KCTsO 

RFAO(S*SOOO) TiTUt 
SOOO FORMAT (lOAA) 

1 CONTINUE 
WRITE (9«90nn) 

WRITE(Af600A) 

0000 FORMAT (aw X) 

9000 FORMAT ClMl///////29X*?4M* * • »0S READ OUT • * *//) 
RFAO(?*2000I NSAMPL*NFXP*STATUS*HATE*NX1 •NX?fNYl*WY?*D 

2000 FORMATf2l2*2Al«An«AA) 

IF(NEXP.E0.99) CALL EXIT 
IFINSAMPL.FO.OI NSamPlsIOO 
WRTTE|9«900l) NEAP«NX1 •NX2«NY1 tNY? 

9001 format fUX*20HEXPEPlMFNT NUMBER s *12/ 

XllXtTMTlME * •2n«lHt«2Il/) 

IF(STATUS.FQ,ATEST(n .0R«STATUS»E0.ATFSTf3») WRITE (9#9002) 
IF(STaTUS.E0.ATESTc2) •0R,STaTUS.EO,ATEST(4) ) write C 9f 9003) 
IF ISTaTUS.FO.ATEST 13) .OR«ST ATUS.EQ.ATFST C*5 ) WRITE <9*9004) 
IF(STaTUS.FO*aTFST(I) , 0R,STATUS,EQ,ATFST<?) ) WRITE <9*9005) 
900? FORMAT <10X*?7H HADIOMFTfcR IN OPERATE MODE) 

9003 format < I OX ♦POM radiometer IN CALIBRATE MODE) 

9004 FORMAT <lOKt?6M TLMPERaTUHFS NOT RFCORDEO/) 

9005 FORMAT < 10 X ♦PPM TEMPERATURES RECOROFD/) 

IF<RATE*E0.ATFST<5>) IRATFsl 
IF<RATF*FO,ATEST<f»)) IRATE*? 

IF<RATE.FQ.ATEST<7) ) IRaTF*4 
1F<RATF.F0.ATEST<»>) IRaTEsS 

WRITE <9*9004) IRATE vNSAMPL 
9000 FORMAT <10X%15M SAMPLE RATE * *n*4M/SFC/ 

X10X*31H NUmRFR of radiometer samples s *13//) 
IF<STaTUS*FQ.ATFST|3).0R,STATUS.E0,ATEST<a)) RO to 30 
READ<?*2001) <IV0LTI1)*I*1*I6) 
pool F0RMAT<l6I4> 

1F(IV0LT<1),NE,9999) 60 TO 10 
WRITEl9*9u07) 

9007 FORMAT <////! OX* 34H CAUrRaTION SEQUENCE WAS NOT FOUND*) 
KCT*0 
60 TO I 
10 CONTINUE 

WRITE <9*900R) < IVOl.T < I) ♦ 1*14* 14) 

900A FORMAT <??X«1SHPDS CAL IHRAT ION/ 11 X ♦ 

X35MZCRO PROBE SCALE FULL SCALE/ 

X11X*I3*9X,T4*1?X*I4//> 


RAOOOOlO 

RAOOOQ20 

RAD00030 

RA000040 

RAOOOOSO 

RA000060 

RA000070 

RAOOOOSO 

RAD00090 

RAOOOlOO 

RAOOOllO 

RAD00120 

RA000130 

RA000140 

RA000150 

RAD00I60 

RA000170 

RAOOOISO 

RAO00I90 

RAD00200 

RAD00210 

RA000220 

RAD00230 

RA00024C 

RA000250 

RAD00260 

RA000270 

RAD00280 

RA000290 

RAD00300 

RAD00310 

RA000320 

RA000330 

RAD00340 

RA000350 

RA000360 

RA000370 

RAD00380 

RAD00390 

RAD00400 

RA000410 

RAD00420 

RAD00430 

RAD00440 

RA000450 

RAD00460 

RAD00470 

RAO004R0 

RA000490 

RA00.500 

RA000510 

RA000520 

RADQ0530 

RA000540 

RA0005SO 
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raOiom fortran pi 


IDSPTI 


OH/U/73 


l?,37,U 


PA6F ^ 


IFdVOLTU**) ,f O.O» IVOtTdSl*! 

no ao tsl«)3 

RES« I VOL T ( n •*»Hb I . / 1 VOL T <1 *> ) 

OFFSFTaO, 

IF<I«F0.1) OFFSET*-. 3 

CALL PRTENPfRtStOFFSET.TEHPRl I» ) 

IF<TFMP«m ,LT.0.» nNPR(I)sRFS«.l 
20 CCSilNUF 

WRITE (4*90091 ( (T ITLE { I • J) • 1*1 * J 0) •TEmPR(j) »J*1*]3) 

9009 FOR AT(2AX«?7H PRORE TFNRERATURFS C0E6 C»// 

X*10X*6H R1 s *10A4*3H * «F6.1.?H C/ 

X*10X*AH R2 s *I0AA*3H s *P#>.1*2H C/ 

X*10X*6H R3 « *I0A«*3H * «F6.1*2H C/ 

X*10X*6H R9 s *10A%*3H s *F6.1*2H C/ 

X«10X«6H R!> * »10A4*3H s *F6.1«?H C/ 

X*10X*OH R6 * *10A4*3H * *Ff>.l*2H C/ 

X*10X*'.H R7 * *10A**3H * *F4.1*2H C/ 

X*10X*AH RR * *10A4«3H * «F6.1*2H C/ 

X«10X*6M R9 s *10A4«3R * *F6.!«?H C/ 

X*10X*7H RlO s *10A4*iH • «F4.1*?H C/ 

X*10X»7H All « «10A4*3H * «F4.1»2H C/ 

X*10X«7M R12 « *lOA4*iH s «F4.1«?H C/ 

X*10X*7H A13 * *10A4*3H s *F4.1,2H C) 

WRITE (3*30AO) NXl*NX2*NYl*NV?*TeMPR 
3000 FORMAT (411 *13F5.1) 

30 CONTINUE 

REAO(2«200?) NUMRAO 
2002 FORMAT (13) 

IF(NUMRA0.EQ.999) go to 39 
READ(2*200n ( f VOLT ( I ) * I«l *NCMRAO) 
if(nomrao.fo.nsampl) go To 40 

39 continue 
WRITE(9.90I0) 

9010 F0»MAT(///10X*24MEWR0P IN RAOlOMETEM OATA) 

KCT»0 

GO TO I 

40 continue 

CALL RAORMSC IvnLT*NUM»AO*XMEAN.HMS*CHI) 

WRITE(9*90ll) XMLANtRMSfCHl 

9011 FORMAT (///10X*2bHMEAN RaUIOMETFR OUTPUT * *F6,1/ 
X10X*16HRMS OEVTATION * «F6.1/ 

X10X«1RhCHI-SOUARE test s «E10.3) 

ICAL*0 

(FISTaTUS.FQ.AUST (2) .0R.STATUS.E».ATEST(4) ) TCAL*! 
WRITE (4*4000) NXl*WX2*NYl*Ny2* ICAL « XMh AN«RMS«CHI 
AnO" FORMAKSIl .?F7.2*t)0.3) 
aCT*KCT*1 

IF(ICAL.EQ.I) go To 100 
KCT*0 

0FN0M*CALR(?) C*« H(1) 

IF (OENnM,LT,.OnOOO: ) OF NOMs, 000001 
APTFMP»53,0S*(XMEAN-rALH( 1) ) /DENOM*TEMPP ( 2) .273,2 
RMSATPsb3,OS*HMS/OFNOK 
WRITE (9*9012) APTEmP,rMvJATP*1TIME 
901? FORMAT (//10X*21MAPAREnT VFMPERATUffE a*F7,2.4M OEG K/ 


AAOOOSAO 

RAOOOS70 

RAO00S80 

RAD00S90 

RAOOOOOO 

RA000610 

RA000620 

RA000630 

RA000640 

HA000050 

RA000660 

RA000670 

RAOOOAOO 

RA000690 

RA000700 

RA000710 

RA000720 

RA000730 

RA000740 

RA0007SO 

RA000760 

RA000770 

RA0007HO 

RA000790 

RAOOOeOO 

RA000810 

RA000820 

RAOOOB30 

RA000840 

RAD00850 

RA000860 

RA000870 

RA000880 

RAD00890 

RAD00900 

RA000910 

RA000920 

RAD00930 

RA000940 

RA0009SO 

RA000960 

RA000970 

RAD00980 

RAD00990 

RADOIOOO 

RAOOlOlO 

RA001020 

RA00103C 

RA001040 

RAOOIOSO 

RAOOIOGO 

RA001070 

RADOIOAO 

RA001090 

RADO) 100 
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RAOIOM FORTRAN PI 


lOsPTl 


0H/l%/73 


U,37.11 
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XlOX*llH«MS noise 

OEG 

«/ 

BAOOlllO 

KlOXtSlHHASFD ON 

CAl IHRATION AT 

TINE s •?Il«lHt«?Il) 

RA001120 

GO TO 1 



NAD01130 

100 CONTINUE 



OA001140 

IF<KCT,FO.n GO 

To 110 


NAOOllSO 

CAL»U>sCALSAV 



RA001160 

calh(?)**rfan 



RAOOllTO 

ITIHE <1>*NX1 



MA001180 

ITlN€<?>sNX? 



RA001190 

ITINE(3J*NV1 



OA001200 

ITINE(A)sNV? 



HA001210 

KCTsO 



RA001220 

GO TO 1 



RA001230 

110 CALSAVsANEAN 



RA0012A0 

GO TO 1 



HA0012SO 

FNO 



RAO012G0 
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I-BTEMP FOPTRAM PI 
VP/CSS 


IbsPTi OH/U/73 

national css* INC. CSTamfORD data CENTER) 


SUPROUT INF PHTE NR I PF S ♦ OF F SF T * Tt PP ) 
niNENSTON ACOF *«C0F IA) *CrOF (6) 
data ACOF/0,Ui?fi^E-OS*0.?OS'»OE-OS,0.3?3<»AC-05* 
X0,AR2?PF-0S*0. f8?30E-0S,0. 1 ISA7F-0*/ 

DATA BC0F/-0,1AA87F-01*»0. 1N33AF -01 »*0.?ASB4E-0 1 • 
X-0.3037ie-0l*-O,iN3l?F-0l*-O,4843AE-0l/ 

DATA CCOF/SI .AOSfSH. 1 10,b3.N4N,A9.?34. 7S.404*B1 .4AS/ 
TCOFsA 

IFtRFS•LT,9*»4.^) CO TO 10 
lFfRFS.6T,3S3P) (iO TO 10 
IFCRES.OT.IPOO) ICOFaS 
lFfRES.GT.1471) IC.>Fs4 
IFfRES.GT.lAlb) 1COFS3 
IFfRES.GT.?P53) ICOFs? 
lFfRES.GT.pAlA) ICOFal 

TFNPaACOF f ICOF ) *RF SR«ES*MCOF f ICOF ) *RES^CCOF f ICOF ) *OFFSFT 

return 

10 TFNPa-l. 

RFTURN 

ENO 


RAGE 1 


RBTOOOlO 
RBT00020 
RBT00030 
RBT00040 
PBTOOOSO 
PBTOOOGO 
PBT00070 
PBTOOOSO 
PBTOOOSO 
PBTOOlOO 
PBTOOllO 
PBTQ0120 
PBT00I30 
P6T00140 
PBTOOlbO 
PBT )01G0 
PBTOOITO 
PBTOOIBO 
PBTOOlOO 
PBT00200 
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PADRmS F0»T»A»I pi TO*MTl OB/U/73 1?,37«S1 PA6C 1 

VP/CSS — national css* INC. CSTANFOPO DATA CENTER} 


SURROUT INF PAOumS(kIN*N*XMEAN«XPMS*ChI) 
OfNENSTON KTNUOO} 

totaled • 

PO 10 7sl*N 
TOTaL«TOTAI *KINai 
10 CONTINUE 

XMEANcTOTAL/N 

total«o. 

DO ?0 1*1 *N 

TOTAL*TOTAL* l*»*EAN-KlNn} 

20 continue 

X*«S3T0TAL/(N<-1 I 
JsO 

TOTALsO. 

00 30 TsUn 
OEV* ( XME AN-K IN ( I > > ••i 
1FC0EV.6T,Q,«XmS} go to 30 
J«J I 

KIN(J)sKlN(T) 

TOTALxTOTAL«KIN(J) 

30 CONTINUE 

xmean«total/j 

T0TAL*0* 

00 40 I*1*J 

TOTAL*TOTAL*(XmEAN-KIn<I» )**2 

40 CONTINUE 

XPNSsSQPI ctotal/*u-im 

LUO 

L?*0 

L3*0 

L4*0 

on 50 I«1*J 

IE(KlNa).6F*(XNEAN*XPMS/P*} I GO TO 41 
IFCKlNd) .6F*XMEAN) 60 TO 42 
IF(KINII) .6F.<XMEAN-XPMS/?*I 160 TO 43 
Ll«LUl 
60 TO 50 

41 L4«L4*I 
60 TO 50 

42 L3*l3*1 
60 TO 50 

43 L2«L?*1 
50 CONTINUE 

01*.30B5*U 

02*.1915»J 

CH1»( CLl-01 ) •42/01 1 ♦ ( CL2«(12) **2/021 
X^ I (L3-021 **2/021 ♦ ( CL4-01 1 **2/0l I 

return 

END 


HAOOOOlO 

RA000020 

RA000030 

RA000040 

PA000050 

RAD00060 

RA000070 

PAOOOOIIO 

RAD00090 

RAOOOlOO 

RAOOOllO 

RA000120 

RA000130 

RA000140 

RA000150 

RA000160 

RA000170 

RA000180 

RA000190 

RA000200 

RA000210 

RA000220 

RA000230 

HA000240 

RA000250 

RA000260 

HAD00270 

RA000280 

RA000290 

RA000300 

RA000310 

RA000320 

RA000330 

RAD00340 

RA000350 

RA000360 

RA000370 

RAD00380 

RAD00390 

RA000400 

HA000410 

RA000420 

RAD00430 

RA000440 

RAD00450 

RA000460 

RA000470 

RAD004*' 

RAOr 



